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i—t \ We calculate the Isgur-Wise function (£q) by measuring the heavy- heavy meson 

transition matrix element on the lattice in the quenched approximation. The standard 
Wilson action is used for both the heavy and the light quarks. Our numerical results 
are compared with various model calculations and experimental data. In particular, 
using a linear fit, we find p 2 = 1.24 ± .26 ± .26. (The slope of £ at the zero-recoil point 
^ is — p 2 .) Using instead the parametrization of £o suggested by Neubert and Rieckert 

and fitting the lattice data up to v ■ v 1 ~ 1.2, we find p^i? = 1.41 db .19 db .19 



Abstract 



1 Introduction 



The physics of heavy quark systems has been one of the active research topics in the last 
few years (for reviews, see Consider a meson system consisting of a heavy quark Q and 
a light quark q. When ttiq — > oo the heavy quark essentially becomes a static color source. 
It can be described by a simple heavy quark effective field theory (HQEFT) with flavor and 
spin symmetries. Because of these new flavor and spin symmetries in HQEFT, the meson 
decay matrix element can be simplified and different decay processes can be related to each 
other. For example, in the case of B — > D decay there are two independent form factors 
in the full theory. However, if we assume both B and D mesons are heavy enough, using 
HQEFT the number of unknown form factor can be reduced to one. We have H 



< D(v')\cy u b\B(v) >= y/m B m D C cb (n)^ (v' ■ v; fi)(v + v' 



(1) 



where v u ,v' u are the four-velocity and mB,rriD are the B and D meson mass, respectively. 
The constant C c b(n) comes from integrating the full QCD contribution from the heavy quark 
mass scale down to a renormalization scale // <C 
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where a(v ■ v') is a slowly varying function of v' ■ v which vanishes at v — v' 0. The Isgur- 
Wise function £q(v' ■ v; jj) represents the interactions between the light degrees of freedom in 
the heavy meson system. It thus has to be calculated using nonperturbative methods. 

On the lattice the heavy meson system can be studied in two different approaches. One 
is to keep the heavy quark dynamical by using the standard Wilson action. This may require 
extrapolation to the physical heavy meson mass of interest. The alternative is to integrate 
out the heavy quark first and derive an effective action including only the light degrees of 
freedom and then perform numerical simulation using this effective action || . Here we stay 
with the first approach (see || for reviews of the second approach). We reported preliminary 
results from our approach at the Lattice '92 conference J|] . 

It is important to note that using the flavor symmetry of HQEFT the Isgur-Wise function 
relevant to the B — ► D decay of Eq. flU) can be obtained also from the D — » D elastic 
scattering matrix element pj 



where 



< D v ,\cq v c\D v >= m D C cc (fi)£ (v ■ v'; fi)(v + v') l 

a(v-v') 
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(3) 



(4) 



This of course requires that the D meson be sufficiently heavy for the onset of the heavy 
quark limit (HQL). 

Conventionally B — » D (here we use B and D as generic names for heavy pseudoscalar 
mesons, they do not necessarily represent the physical B and D mesons) transition matrix 
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elements can be parametrized as 

< D(p')\v v \B(p) >= U(g 2 W + p)u + f-(q 2 )(p-p') 



(5) 



where q 2 = (p' —p) 2 is the momentum transfer between the initial and final states and V v is 
a vector current. One can also use an equivalent form 



< D(p')\v u \B(p) >=f + (q 2 
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with relation between the form factors /+, /_ and fo 

(ml - m 2 D )f (q 2 ) = (m| - m 2 D )f + (q 2 ) + q 2 f^{q 2 ) . 



(6) 



(7) 



For the elastic scattering process (tub = ttid), we have f-(q 2 ) = 0, which follows from 
current conservation. Using the relation v v = p u /m, Eq. (13) becomes 



< D(y')\V„\D(y) >= m D f + (q 2 )(v' + v) u . 
Comparing this with Eq. (||) one finds the simple relation between /+ and 



(8) 



(9) 



The lattice calculation method for /+ has been well established We will use that 

method through equation (9) and deduce £o- 



2 Considerations in the Lattice Calculation 



Choosing the lattice parameters. HQEFT becomes valid when the heavy quark mass mg ^> 
Aqcd- Therefore we need to choose the lattice parameters to satisfy this constraint. For 
example, at (3 — 6.0 the inverse lattice spacing a -1 w 2.0 GeV. If we take the physical value 
Aqqjj ~ 0.2 GeV, then in lattice units Aqcd = 0.1. We need to choose the heavy quark mass 
amQ to be much larger than 0.1. However, one would expect large lattice spacing artifacts 
for arriQ>l. At (3 = 6.0 this means that the physical mass of the heavy quark is limited to 
be in the range of the D meson or smaller. At larger (3 values, the lattice spacing is smaller, 
and we can therefore accommodate mesons with larger physical masses on the lattice. In 
this work we choose the heavy quark mass in the range 1.5 — 3 GeV by tuning the hopping 
parameter kq for the heavy quark Q. On the other hand, the light quark q in the heavy 
meson should have a mass m q < Aqcd- This can be achieved by changing the hopping 
parameter n q for the light quark and then extrapolating to the chiral limit K q)Cr , where the 
extrapolated pion mass is approximately matched to the physical pion mass m n = 140 MeV. 

How far can v-v' change on the lattice? In the lattice calculations for the matrix element, 
we always have either the initial or the final particle at rest. Thus 

„ ? /_ J E D /m D if u' = (0,0,0, 1) , . 

1 E' D /m D if v = (0,0, 0,1) ' [W) 



2 



with Ed = \ m D + P 2 - F° r a spatial lattice L = 24, we have injected momenta 

p=— (1,0,0) ,—(1,1,0) . (11) 

Since mo « 1.0 in lattice units, one gets t> • t>' ~ 1.034 and 1.066 respectively. To get larger 
values for v ■ v' one needs to inject larger lattice momenta which would in turn introduce 
large statistical noise in the matrix element calculations. Thus v ■ v' can not be much larger 
than ~ 1.1 in current simulations. 

Removing the lattice artifacts; normalization. Since the heavy meson mass is near 1 in 
lattice units the lattice artifacts could be significant. Ultimately the lattice artifacts can be 
brought under control either by comparing data at different (3 values or by using improved 
lattice actions. However, for simulations at a given /3 value there are several ways to check 
the size of the lattice artifacts. Also, for the purpose of this work, we are able to get rid of 
most of the undesirable effects by imposing the conservation of heavy quark flavor current, 
as we will explain below. 

For the elastic scattering process the matrix element has the simple form given in Eq. 
(^) . This is derived from Eq. (^j) which is based on the assumption that the theory is Lorentz 
invariant. The lattice theory, however, is not exactly Euclidean rotational invariant due to 
the finite lattice spacing a. Thus /_ is not exactly zero. The amplitude of /_ (or /_//+) 
gives a measure for the violation of the Euclidean invariance on the lattice. 

We can also estimate the size of the lattice artifacts by checking the simulation results 
against known continuum matrix element values at some special points. For example, when 
both the initial and the final D mesons are at rest, v = v' = (0, 0, 0, 1), the continuum matrix 
elements of C74C is known because of the quark flavor current conservation 

< D\c^c\D >= 2m D . (12) 

At this so-called "zero-recoil point" we have 

£o(l) = 1 • (13) 



Both Eqs. (^j) and ( 12 ) will have O(a) corrections on the lattice which can come from different 
origins. 

Part of the 0(a) effect can be disposed of, to a good approximation, by including a 
normalization factor 

< VW(0) >cont= 2Ku e m < ^(x)^(0) > latt , (14) 



where 



with u the "tadpole improvement" factor. 



e m = 1 + — f ^- - ) , (15) 
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Another correction comes from the use of (nonconserved) local vector current V v = c^ v c 
on the lattice ||. This effect can be corrected by introducing a rescaling factor Zy c in the 
vector current. In perturbation theory Z l y c is calculated to be || 

Z v = 1 " 27-5^. (16) 



In general the size of lattice artifacts will be momentum dependent so the 0(a) correction 
will be different for Eq. (^) and Eq. (§). However, since we are using only small momentum 
injections, we assume that the leading 0(a) correction is approximately independent of p,p' 
and the Lorentz index v. Thus by enforcing Eq. ( |12|) as the normalization condition for 
the matrix element < D(v')\V v \D(v) >, corrections due to both the exp(m) factors and Zy c 
factors are taken care of automatically (together with any other multiplicative 0(a) factors). 

Corrections due to finite heavy quark mass. At what mass scale the heavy quark limit sets 
in is an important question. The charm quark mass is not that much larger than Aqcd- If we 
choose the p- meson mass m p as the typical scale for QCD, it becomes even less clear if m c is 
heavy enough. In lattice calculations of the matrix elements the size of the corrections to the 
symmetry limit (i.e., itlq — > oo) have been observed to be process dependent. For example, 
the heavy meson decay constant fp is found to have a mass correction that is as large as 
50% at .D-meson mass scale fTO] , |TTJ ; while the form factor Ai(q 2 ) at the zero- recoil point 
q 2 = q 2 max for D — > K* decay shows very little mass dependence even for fairly small values 
of the heavy quark mass [TJj . Both observations are consistent with theoretical expectations. 
In HQEFT the leading correction due to the heavy quark mass is in general 0(1/ 'mg), as is 
the case for fp. However, in the case of D — > K*, the 0(1/ mo) correction vanishes exactly 
at the zero-recoil point and the finite mass correction there is at most 0(1 /niq) |13|. Thus, in 
lattice calculations, the observations of the linear 0(1/7tiq) correction for fp calculation and 
little mass dependence of A^q^^) for D — > K* decay seem to indicate that the finite mass 
corrections are under control within HQEFT at a scale ~ m c and that numerical simulation 
results in this range can be used and extrapolated reliably to obtain the heavy quark limit 
results. 

In principle, for data obtained at fixed masses one may estimate the size of 0(1 /ttiq) 
corrections in the following way: the matrix element can be in general written as 

< A(v')\V u \B(v) >= ^^Cab [f + (v ■ v')(v' + v) v + /> ■ v')(v' - v) v ] . (17) 

This parametrization is equivalent to the conventional parametrization Eq. (|5|); we have 
relations 

1 lm. a Im.n 

(18) 
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f- 




m B V m A 



m B 



m A 



f- 



(19) 



In the m<3 



oo limit because of the constraints on the effective heavy quark fields 

luVuK = h v , h v j u v u = h v , (20) 
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we have [|TJ 

/_ = . (21) 
At finite tuq, /_ is nonvanishing. Thus the value of /_ or the ratio 

f_ = 1 + /-//+ -%*(!-/-//+) 

/+ 1 + /-//+ + -/-//+) ' 1 ] 

gives a measure of the 0(1 /mo) correction to the leading order HQEFT. Unfortunately, for 
the elastic scattering process considered in this paper this relation is not useful because in 
this case we have itla = mg so that /_//+ = /-//+ = in the continuum. However, in 
principle, equation fl2"2p may be used to check the finite mass correction in computations of 
B — > D decay. 



3 Numerical results 



We use the standard Wilson action for the propagation of both heavy and light quarks and 
work in the quenched approximation. The gauge configurations used are listed in Table 1. 
The techniques for measuring the two-point function and the three-point matrix elements 
are standard H [J. Here we only list the results for the measured uid, /+(0) and f+(q 2 ) 
in Tables 2-7. To increase the statistics, we have used symmetry properties of the Green 
functions and averaged over ±t and ±p whenever possible. 

Estimating the size of the lattice artifacts. Data for /+(0) and f+(q 2 ) listed in Tables 2-7 
are direct numerical results. They include all the lattice artifacts. For example, comparing 
to the known continuum value /+(0) = 1, we observe that the lattice artifacts at /3 — 6.0 are 
typically 20% — 40% at kq = 0.118 and less than 10% at k,q = 0.135. From the computation 
of the ratio /_//+ we find the violation of Euclidean invariance is typically 5% — 15% for 
kq = 0.118 and 3% — 10% for k,q = 0.135. Note that the reduction of the lattice artifacts 
when kq is changed from 0.118 to 0.135 agrees with our intuitive expectations that those 
lattice artifacts that are due to the heavy quark mass approaching the lattice ultra-violet 
cutoff should decrease with a decrease in mass of the heavy quark. One may try to use the 
factors Zy c and e m to remove some of the artifacts. At /3 — 6.0 we get from the perturbative 
calculation, Eq. (|T|), that Z l J c ^0.7 (using the "boosted" effective gauge coupling g 2 « 1.7 
as suggested in ref flip . The factor e m is about 2 for kq = .118 and 1.5 for kq = .135. 
Including both these factors the corrected /+(0) become around one within errors. This is 
in agreement with other observations |IJ that Zy c and e m factors seem to account for the 
largest part of the lattice artifacts. 

Normalization and the Isgur-Wise function. As explained in Section 2, we define ^q(v • 
v') = f + (q 2 )/f + (0) in an effort to correct for possible lattice artifacts. We list the results 
in the last two columns of Tables 2-7 []. This method for deducing £ has the significant 

1 We have also computations at (3 — 5.7 with results basically consistent with data listed here within large 
errors. However, the heavy quark mass used (corresponding to k — 0.094) is extremely high for that j3 with 
potentially very large lattice artifacts. They are therefore not used in this work. 
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advantage that it is free of uncertainties in factors such as Zy c and e m . Note that in Eq. 
(|9]) there is a factor C cc in the connection between / + and £q. This factor comes from 
integrating out the QCD effects from the heavy quark scale down to a light scale \x. For the 
current calculation, however, the heavy quark and hard gluon effects in QCD are included 
dynamically, and we are therefore determining the combination C cc £o at arbitrary scale /x. 
For simplicity we choose fi = mo and set C cc = 1 according to Eq. (f|). 

Finite volume effects and the dependence on the light quark mass. In an effort to keep 
finite volume effects under control, we first study our data with light quark mass (m q ) in 
the range of the s-quark mass (m s ). On the lattice we choose k = .154 at (3 — 6.0 and 
k = 0.149 at (3 = 6.3. The physical meson masses at these (3 and k values are in the range 
of 0.5-1 GeV, corresponding to bound states of ss. The results for the Isgur-Wise function 
are plotted in Fig. 1. All data seem to fall on a smooth curve. In particular, data for the 16 3 
and 24 3 lattices at (3 = 6.0 agree to a very good approximation (i.e., well within statistical 
errors ) assuring us that in this limit finite size effects on these lattices are small. Since the 
physical volume of the 24 3 lattice at (3 = 6.3 is the same as the physical volume of the 16 3 
lattice at f3 = 6.0, we expect the finite size effects on the 6.3 lattice also to be small. 

As discussed in Section 2, for the elastic scattering of physical D mesons, one should 
extrapolate the light quark hopping parameter to the chiral limit k — *> n cr . By taking this 
limit on a finite lattice, we expect not only a shift of the Isgur-Wise function due to the 
change of the light quark mass, but also a finite volume effect. The former effect is physical 
but the latter one is due to the limitations of a finite lattice size and should be brought under 
control by carefully comparing results on different lattice sizes. In practice, however, it is 
hard with a limited set of data to separate the light quark mass dependence from the finite 
size effect. Since, as Table 8 shows, the physical linear size of all our lattices are roughly 
in the range of (100 MeV)" 1 , as the light quark mass gets smaller than m s we expect finite 
size effects to become important, especially on the two smaller lattices. Indeed, comparison 
of Fig. 2 with Fig. 1 shows that all the data is shifted up in the chiral limit. The size 
of shift is the smallest on the (3 = 6.0, 24 3 lattice, which has the largest physical volume. 
Therefore, the (3 = 6.0, 24 3 lattice is the most suitable one for extracting the dependence of 
the Isgur-Wise function on the light quark mass. Comparison of the numbers in Tables 4 
and 5 corresponding to k = 0.154 with those from the extrapolated case (i.e., for k = 0.157) 
indicates that the form factors and the ratio / + // + (0) change very little (i.e., < 2%) at a 
fixed momentum p. Most of the change in the Isgur-Wise function comes from the change 
of the D meson mass (with the light quark mass) which in turn shifts the value of v ■ v' . We 
note that the light quark mass dependence of the Isgur-Wise function has also been studied 
using chiral perturbation theory |L5] and model calculations |16[ . It was found that the shift 
due to the light quark mass is small: at y — 1.2 the shift in £o ranges from « —2% to about 
7%. 

Comparison with model calculations and experimental data. As discussed in the previous 
two paragraphs, we quote the m q ~ m s data as our final results. They are listed in Tables 
2-7 (At (3 = 6.0, n q = 0.154 on the 16 3 and 24 3 lattices and at (3 = 6.3, n q = 0.149 on the 24 3 
lattice) and plotted in Fig. 1. For comparison we also plotted the theoretical upper bound 
(dashed line) on the Isgur-Wise function derived from algebraic sum rules by Bjorken |~7 . 
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The lower curve (dotted line) is the lower boundQ derived in Ref . | IS | . All of our data appear 
to be, within errors, below the upper bound of Bjorken. The solid curve going through the 
data points is a fit to the form 



-(2p^-i; 



v' . (23) 



y + 1 

Taking into account the correlations in the data, we find, after fitting, p 2 NR = 1.41(19) with 
X 2 /dof ~ 13/12. A similar fit to the data (which, as explained in the preceding paragraph, 
is in this case only from the (3 = 6.0, 24 3 lattice) in the chiral limit gives p 2 NR = 1.09(28) with 
X 2 /dof ~ 1.4/4. However, this shift between the central values (1.41 vs. 1.09) cannot be 
taken as a reliable reflection of the light quark mass dependence of the Isgur-Wise function. 
Indeed, holding p 2 NR = 1.41 fixed also provides a good fit ( with \ 2 ~ 2.6/4 ) to the 
(3 = 6.0, 24 3 data in the chiral limit. Thus it is reasonable to assume that the shift of the 
Isgur-Wise function from m q = m s to the chiral limit is within the quoted statistical errors. 
So far as other systematic effects such as the residual 0(a) effects and 0(l/rriQ) effects are 
concerned, we expect them to be rather small, since we are able to get a good universal fit 
to all the data. We therefore estimate the total systematic error to be at most the same size 
as the statistical error, and quote p 2 NR = 1.41 ± .19 ± .19 as our result for m q ~ m s as well 
as in the chiral limit. 

Close to v ■ v' — 1 the Isgur-Wise function can be parametrized as 

£o(v V) = 1 - p 2 (v ■ v 1 - 1) + 0((v ■ v' - l) 2 ) . (24) 

If we fit our data points to this linear form for v ■ v ' < 1.06 we get p 2 = 1.24(26) with 
X 2 /dof = 2.6/5. Once again, that value of p 2 also provides an acceptable fit (x 2 ~ 1.8/3) 
for the three data points with v ■ v ' < 1.06 (of the (3 = 6.0, 24 3 lattice) in the chiral limit. 
Therefore for the linear fit we quote the result as p 2 = 1.24 ± .26 ± .26. We list our linear 
fit result in Table 9 along with p 2 values estimated by other authors. We note that lattice 
and various model calculations seem to be in rough agreement within errors. However, the 
value of Blok and Shifman tends to be rather low compared to the lattice result. 

In Fig. 3 we plot our lattice results taken from Fig. 1 against experimental data from AR- 



GUS [p6|l . We have normalized the experimental graph with respect to \ V c b\y t b /(1.32 psec.) = 
0.047. For comparison we also plotted the theoretical upper bound |I7[ (dashed line), our 
fit to Eq. (p3|) (solid line) and the curve obtained in Ref. [18| (dotted line). It is gratifying 



that the lattice data has much smaller errors than the experimental results in the region 
v ■ v' < 1.15. 



4 Conclusions 



We have computed the Isgur-Wise function on the lattice from the elastic scattering matrix 
element. The calculation is done using the standard Wilson fermion formulation for both 

2 Recent work brings into question the validity of this lower bound for the Isgur-Wise function p9| . 
However, one may take the result of Ref. [pj| as a model calculation. 
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heavy and light quarks in the quenched approximation. We find that all of our numerical 
results are within the theoretical upper bound JHJ. The Isgur-Wise function (£ ) is deduced 



by taking ratios of form factors. In this approach the lattice artifacts appear to be under 
control. Although the heavy meson mass in our calculation is typically order one in lattice 
units, the residual 0(a) effects are apparently small, as indicated by direct comparison of 
the data at /3 — 6.0 and (3 = 6.3. Further checks can also be done with an improved action 
scheme ^7f . The lattice meson mass we used is in the range of physical D meson. Thus there 



are potential 0(1 /mo) corrections. Within the limited range accessible in our simulations 
we have not observed strong rriQ dependence. 

Our best lattice results are obtained when the light quark mass is set at the scale of 
s-quark mass. Data from our largest lattice seems to indicate that the Isgur-Wise function 
changes by less than our statistical error between m q ~ m s and the chiral limit. A model 
independent linear fit using data close to v ■ v' = 1 gives the slope p 2 = 1.24 ± .26 ± .26. 
Comparison of the slope p 2 with continuum model calculations are given in Table 9. A direct 
comparison to the ARGUS experimental results |26j is given in Fig. 3. Our results are also 
in good agreement with the preliminary results of UKQCD Collaboration P7[ . 
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p 


L 3 T 


# conf. 


fitting range 


6.0 


16 3 39 


19 


10 < t < 15 


6.0 


24 3 39 


8 


10 < t < 15 


6.3 


24 3 61 


20 


15 < t < 25 



Table 1: List of gauge configurations used in this work. 



Kg 




/+(0) 


V 


/+ 


v ■ v' 


£o 


0.152 


1.237(6) 


0.69(12) 


(1,0,0) 
(1,1,0) 


0.65(12) 
0.62(14) 


1.0408(4) 
1.0790(8) 


0.953(22) 
0.911(57) 


0.154 


1.202(8) 


0.74(15) 


(1,0,0) 
(1,1,0) 


0.72(16) 
0.70(18) 


1.0434(5) 
1.0842(10) 


0.974(30) 
0.945(79) 


0.155 


1.186(8) 


0.77(19) 


(1,0,0) 
(1,1,0) 


0.76(21) 
0.74(24) 


1.0445(6) 
1.0869(11) 


0.983(42) 
0.95(11) 


K cr 


1.151(8) 


0.82(19) 


(1,0,0) 
(1,1,0) 


0.82(21) 
0.81(24) 


1.0567(8) 
1.1103(15) 


1.00(4) 
0.99(11) 



Table 2: Data and extrapolation to k ct on the /3 = 6.0, 16 3 x 39 lattice. The heavy quark 
hopping parameter is set to kq = 0.118. k ct 0.157. 



Kg 


m D 


/ + (0) 


V 


u 


v ■ v' 


^0 


0.152 


0.885(6) 


1.00(12) 


(1,0,0) 
(1,1,0) 


0.89(12) 
0.77(14) 


1.085(1) 
1.160(2) 


0.883(29) 
0.771(72) 


0.154 


0.847(7) 


1.03(15) 


(1,0,0) 
(1,1,0) 


0.92(16) 
0.79(18) 


1.0937(14) 
1.1749(25) 


0.894(40) 
0.77(9) 


0.155 


0.829(7) 


1.06(19) 


(1,0,0) 
(1,1,0) 


0.94(20) 
0.78(22) 


1.0980(15) 
1.1841(27) 


0.894(40) 
0.74(12) 


K cr 


0.791(7) 


1.08(19) 


(1,0,0) 
(1,1,0) 


0.98(20) 
0.79(22) 


1.116(2) 
1.222(4) 


0.90(5) 
0.73(12) 



Table 3: Data and extrapolation to k ct on the (3 = 6.0, 16 3 x 39 lattice. The heavy quark 
hopping parameter is set to kq = 0.135. k ct ps 0.157. 
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Kq 


m D 


/+(0) 


V 


/+ 


v ■ v' 




0.152 


1.241(8) 


0.61(13) 


(1,0,0) 
(1,1,0) 


0.60(14) 
0.57(14) 


1.0185(2) 
1.0366(4) 


0.973(20) 
0.934(35) 


0.154 


1.202(10) 


0.64(13) 


(1,0,0) 
(1,1,0) 


0.63(14) 
0.60(14) 


1.0201(3) 
1.0396(7) 


0.975(20) 
0.938(42) 


0.155 


1.182(11) 


0.66(12) 


(1,0,0) 
(1,1,0) 


0.64(12) 
0.62(12) 


1.0209(4) 
1.0412(8) 


0.973(20) 
0.936(49) 


HL cr 


1.144(10) 


0.69(13) 


(1,0,0) 
(1,1,0) 


0.67(14) 
0.65(14) 


1.0258(4) 
1.0511(9) 


0.98(2) 
0.94(5) 



Table 4: Data and extrapolation to n cr on the (3 = 6.0, 24 3 x 39 lattice. The heavy quark 
hopping parameter is set to kq = 0.118. K cr ~ 0.157. 



K q 


m D 


/ + (0) 


V 


/+ 


v ■ v' 


£o 


0.152 


0.889(7) 


0.96(10) 


(1,0,0) 
(1,1,0) 


0.90(10) 
0.84(10) 


1.0413(6) 
1.0802(12) 


0.944(17) 
0.873(33) 


0.154 


0.847(8) 


0.99(10) 


(1,0,0) 
(1,1,0) 


0.94(10) 
0.87(11) 


1.0458(8) 
1.089(3) 


0.949(19) 
0.883(40) 


0.155 


0.826(9) 


1.01(10) 


(1,0,0) 
(1,1,0) 


0.96(10) 
0.90(11) 


1.0476(10) 
1.093(5) 


0.948(18) 
0.884(43) 




0.784(9) 


1.05(10) 


(1,0,0) 
(1,1,0) 


1.00(10) 
0.90(11) 


1.0543(12) 
1.106(2) 


0.95(2) 
0.89(4) 



Table 5: Data and extrapolation to n cr on the (3 = 6.0, 24 3 x 39 lattice. The heavy quark 
hopping parameter is set to kq = 0.135. K cr ~ 0.157. 
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K q 


m D 


A(o) 


V 


/+ 


v ■ v' 




0.148 


0.928(4) 


0.74(8) 


(1,0,0) 
(1,1,0) 


0.69(8) 
0.66(9) 


1.0350(3) 
1.0708(5) 


0.929(14) 
0.89(4) 


0.149 


0.906(4) 


0.78(10) 


(1,0,0) 
(1,1,0) 


0.73(10) 
0.71(11) 


1.0373(3) 
1.0752(5) 


0.937(16) 
0.91(5) 


0.150 


0.884(4) 


0.86(13) 


(1,0,0) 
(1,1,0) 


0.83(13) 
0.84(17) 


1.0401(3) 
1.0805(7) 


0.960(23) 
0.98(9) 


0.1507 


0.868(9) 


0.99(18) 


(1,0,0) 
(1,1,0) 


0.99(19) 
1.05(27) 


1.0423(6) 
1.0836(17) 


0.993(32) 
1.06(12) 




0.848(9) 


1.02(18) 


(1,0,0) 
(1,1,0) 


1.04(27) 
1.12(27) 


1.0466(6) 
1.0912(12) 


1.02(3) 
1.12(12) 



Table 6: Data and extrapolation to n cr on the = 6.3, 24 3 x 61 lattice. The heavy quark 
hopping parameter is set to kq = 0.125. K cr ~ 0.1516. 



Kq 


m D 


/+(0) 


V 


u 


v ■ v' 


Co 


0.148 


0.587(2) 


1.16(7) 


(1,0,0) 
(1,1,0) 


0.99(6) 
0.86(12) 


1.0894(6) 
1.1743(10) 


0.85(3) 
0.74(11) 


0.149 


0.562(3) 


1.18(8) 


(1,0,0) 
(1,1,0) 


1.05(8) 
0.98(16) 


1.0994(11) 
1.193(2) 


0.89(3) 
0.78(12) 


0.150 


0.537(3) 


1.22(12) 


(1,0,0) 
(1,1,0) 


1.09(11) 
1.06(21) 


1.1128(14) 
1.216(3) 


0.89(5) 
0.87(19) 


0.1507 


0.519(4) 


1.27(18) 


(1,0,0) 
(1,1,0) 


1.20(17) 
1.33(34) 


1.125(2) 
1.235(3) 


0.95(6) 
1.05(28) 




0.496(4) 


1.31(18) 


(1,0,0) 
(1,1,0) 


1.27(17) 
1.48(34) 


1.131(2) 
1.249(4) 


0.98(6) 
1.08(28) 



Table 7: Data and extrapolation to n cr on the = 6.3, 24 3 x 61 lattice. The heavy quark 
hopping parameter is set to kq = 0.140. K cr ~ 0.1516. 
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p 


a" 1 


L 


( a L) p hys 


K Q 


m D 


6.0 


2.0 GeV 


16 


125 MeV 


0.118 


2.3 GeV 


6.0 


2.0 GeV 


24 


83 MeV 


0.118 


2.3 GeV 


6.0 


2.0 GeV 


24 


83 MeV 


0.135 


1.6 GeV 


6.3 


3.2 GeV 


24 


125 MeV 


0.125 


2.7 GeV 


6.3 


3.2 GeV 


24 


125 MeV 


0.140 


1.3 GeV 



Table 8: Physical values of the lattice spacing a, lattice volume aL and "D" meson mass 



Lattice 


1.24 ± .26 ± .26 


Rosner [21] 


1.59(43) 


Neubert @ 


1.42(60) 


Burdman |23 


1.08(10) 


Jin, Huang and Dai |24 


1.05(20) 


Blok and Shifman |25| 


0.65(15) 


Bjorken [0 


0.25 < p 2 


de Rafael and Taron |18 


p 2 < 1.42 



Table 9: Comparison of the lattice calculation, using a linear fit to £o, with various model 
calculations for p 2 at v ■ v' = 1. 
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Figure Captions 



Fig. 1: The Isgur-Wise function is plotted against v ■ v'. The lattice data are obtained for 
m q ~ m s- The parameters are: (3 = 6.0, kq = 0.118, K q = 0.154 on 16 3 x 39 lattice (open 
triangle), (3 = 6.0, kq = 0.135, K q = 0.154 on 16 3 x 39 lattice (solid triangle), (3 = 6.0, kq = 
0.118, K q = 0.154 on 24 3 x 39 lattice (open circle), (3 = 6.0, k q = 0.135, K q = 0.154 on 24 3 x 39 
lattice (solid circle), (3 = 6.3, kq = 0.125, K q = 0.149 on 24 3 x 61 lattice (open square), and 
(3 = 6.3, kq = 0.140, K q = 0.149 on 24 3 x 61 lattice (solid square). Only statistical errors are 
shown. The horizontal errors on the data points are smaller than symbol sizes. The dashed 



curve is the upper bound on the Isgur-Wise function [ 17fl . The solid curve is a fit of the 



lattice data to Eq. fl2"3"|). The dotted curve is derived in Ref. 



Fig. 2: The data for the Isgur-Wise function obtained after extrapolations to the chiral 
limit. The parameters are: (3 = 6.0, kq = 0.118, on 16 3 x 39 lattice (open triangle), (3 = 
6.0, k q = 0.135, on 16 3 x 39 lattice (solid triangle), (3 = 6.0, k q = 0.118, on 24 3 x 39 lattice 
(open circle), (3 = 6.0, kq = 0.135, on 24 3 x 39 lattice (solid circle), (3 = 6.3, kq = 0.125, on 
24 3 x 61 lattice (open square), and (3 = 6.3, kq = 0.140, on 24 3 x 61 lattice (solid square). 
The data points can be compared to the corresponding ones in Fig. 1. For reference we 



have shown the upper bound (dashed curve) on the Isgur-Wise function [17] and the curve 
(dotted) from Ref. 



Fig. 3: Comparison with experimental data. The open circles are ARGUS experiment data 



from Ref. [gfj (We have normalized the graph with respect to \V c b\y tb/1-32ps = 0.047). 
The crosses are the lattice results from Fig. 1. The solid, dashed and dotted curves have 
the same meaning as in Fig. 1. The horizontal errors on the lattice data points are smaller 
than symbol sizes. (Note that the experimental point around v ■ v' — 1.1 is shifted slightly 
for clarity.) 
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